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ABSTRACT 



An analytical model was constructed to predict the 
response of the ball obturator in a ball obturated, spin 
stabilized tubular projectile (BOT) . Torques acting on 
the obturator were considered to be due to fluid and sliding 
friction forces. An experimental apparatus was designed, 
and built to simulate the BOT. Data obtained were the times 
required for the obturator to nutate through a known angle 
over a range of projectile angular velocities. Utilizing 
published coefficient of sliding friction information and a 
linear approximation of the exact solution, the trend in the 
data to go through a minimum value was duplicated. A least 
squares fit of the linear approximation to the data allowed 
the formulation of a new sliding friction coefficient, 

Using this value of Pg , a plot of the least squares fit, 
values from the exact solution and the experimental data 
was constructed. The correlation is extremely good so that 
confidence in the predictive capabilities of the model was 
gained. 
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N0r4ENCLATURE 



(Only primary symbols are listed. Intermediate quantities 
are defined in the text.) 

A major mass moment in inertia, N-m-s 

C minor mass moment of inertia, N-m-s 

G applied torque, N-m 

h width of gap between ball and projectile, m 

M magnitude of applied torque, N-m 

r radius of hole through projectile, m 

rg radius of hole through ball, m 

R radius of ball, m 

6 sin"^ (rs/Rg) 

d,4i,ip Euler angles (see Fig. 7) 

X (C/A)-l 

Vij fluid viscosity (Pa-s) 

tg coefficient of sliding friction, dimensionless 

w angular velocity 

Subscripts 

cr refers to condition when 0=0cr“2s 

o initial value 

p projectile 

r refers to relative angular velocity 

s sphere (ball obturator) 
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x,y,z refers to inertial coordinate system 
1,2,3 refers to body-fixed coordinate system 

Superscripts 

• denotes differential with respect to time 
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I. INTRODUCTION 



A. BACKGROUND 

The tubular projectile is not a new concept. As early as 
1858, the usefulness of tubular projectiles was recognized for 
specific purposes. In that year, Joseph Whitworth (Whitworth 
Threads) wrote about and included an illustration of the tubular 
projectile in the section on Rifled Firearms of his Miscellaneous 
Papers on Mechanical Sx±ijects [1]^. The projectile pictured 
was hexagonal with a circular hole. Whitworth noted its par- 
ticular "...effectiveness in perforating elastic materials 
which prevented them from closing up." The projectile utilized 
a wooden sabot. Whitworth also said that the tubular projectile 
penetrated deeper into masonry than any he was acquainted with. 

The Krnka-Hebler projectile was the next projectile refer- 
enced chronologically [2]. According to Charters and Thomas, 
this projectile was reported in the Allgemeine Schweitzerische 
Militarzeitung as having been so successful that the Ordnance 
Department in the United States carried out firing tests of 
caliber .30 tubular projectiles in 1894. A description of 
those test firings and an updating of the results was performed 
by I.E. Segal and included in Ref. 2. In the 1894 tests, a 
vulcanized fiber sabot was used to push the projectile out the 
barrel . 



Numbers appearing 



in brackets refer to the list of 



re ferences . 
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The purpose of the 1894 test was to compare trajectories 
of standard and tubular projectiles. This was done by com- 
paring the vertical drop on target at a given range. 

Segal's report in Ref. 2 indicates that the drag co- 
efficients computed from the 1894 results agree closely with 
those of Charters and Thomas [2], even though their report 
was published fifty years later (1944) and their experimental 
apparatus was more sophisticated. 

B. RECENT DEVELOPMENTAL STUDIES 

The results obtained by independent researchers as to the 
value of the tubular projectile as an alternative to the con- 
ventional round are contradictory. It seems as if there was 
as much research devoted to discrediting the tubular pro- 
jectile as there was to investigating its merits. Frank and 
McLaughlin [3] have accumulated a great deal of data from 
various sources in an attempt to "objectively compare" the 
merits of the tubular and "conventional shapes". The authors 
concluded that the tubular has no particular advantages over 
"well designed" conventional projectile shapes. 

Their findings disagree with most of those discussed in 
the following sections . 

1. Range Tests 

Winchenbach, Daniel and Edgar [4] conducted range 
tests of six configurations of tubular projectiles and con- 
cluded that the drag coefficients were significantly lower 
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than the standard High Explosive Incendiary (HEI) projectile 
of the same caliber. Only projectiles of the same bore size 
were compared. The experimental models were constructed from 
standard 20-mm ammunition by boring various hole sizes to 
obtain ”... area ratios (A^/Aj^) of 0.7, 0.8, 0.9 and 1.0.” 
Where A-j. is the throat area and Aj^ is the inlet area of the 
projectile. The drag coefficients of the projectiles with 
area ratios of 0.8, 0.9 and 1.0 were less than half that of 
the projectile with area ratio of 0.7. The higher drag co- 
efficient was attributed to the underexpanded flow at the 
exit resulting in high base drag. 

Range tests conducted for the concept evaluation of 
the 20-mm tubular projectile for the Vulcan Weapons System [5] 
showed that the tubular projectile performed better than the 
standard HEI round, with lower drag, shorter time of flight for 
a given range (30 percent shorter at 1000 meters, 40 percent 
shorter at 2000 meters) and similar dispersion characteristics 
on target. 

Recent tests at the Naval Weapons Center, China Lake, 
California of the Ball Obturated Tubular Projectile (BOT) 
have shown a definite advantage in that the low drag of the 
txibular projectile means slower retardation of velocity and 
hence a higher terminal kinetic energy [6]. (The BOT is the 
object of this investigation and will be described in detail 
in later sections of this thesis.) 
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2 . Target Impact Tests 

The higher the energy delivered to the target, the 
greater the damage. Since kinetic energy is directly propor- 
tional to the first power of mass and to the square of the 
velocity, the reduced mass of the tubular projectile is com- 
pensated by a higher muzzle velocity and therefore a potential 
for higher kinetic energy. However, due to the improved drag 
characteristics, the reduction in mass to obtain the higher 
muzzle velocities may not be necessary, depending upon the 
imposed performance criteria. 

Target impact studies have been carried out with 
various projectile configurations and target types. 

Rhethorst, et.al. [7] conducted impact studies of 7.62-itvm 
tubular projectiles on helmets. Tests showed that even with 
the same energy, the tubular projectile penetrated further. 

Kitchen and Keeser [8] conducted studies for the 
Air Force on the impact effectiveness of tubular projectiles 
on simulated aircraft fuel cells. These tests were conducted 
with steel and depleted uraniiam (DU) tubular projectiles and 
standard 20-mm HEI projectiles. The projectiles were fired 
at dox±)le panels at varying degrees of obliquity. Of the 
forty-seven tests of the steel tubular projectiles which im- 
pacted the target, twenty breached the rear panel up to 
angles of 70*^. The DU projectiles breached the rear panel 
even at angles of 85° and fires were started in three of the 
five DU tests. The standard M56 HEI projectiles failed to 
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breach the rear panel in each of the five firings even though 
severe damage to the front panel and two fires occurred. 



Brunsvold and Kalivretenos [9] conducted a program 
to test the effectiveness of the tubular projectile against 
a simulated cruise missile warhead. Tests were conducted 
with 20 and 40-mm rounds. Results of those tests are 
classified, and therefore not presented here. It was noted 
[9] that due to the improved stability characteristics of 
the tubular projectile, only half the spin rate need be im- 
parted by the rifling in the barrel to obtain the same sta- 
bility as a standard projectile, which indicates less barrel 
wear. 

3. Weapons System Capability 

For the tubular projectile to be practical, it must 
be able to interface with existing weapons systems. Refer- 
ence 5 is a report of the 1978 tests conducted by the Army 
to evaluate the tubular projectile in the Vulcan weapons 
system. 

Because the nose of the tubular projectile is flat 
compared to the standard projectile (See Fig. 1), it " . . . 
did not lend itself to chambering in the weapon.” Personnel 
safety required that the weapon be remotely operated and 
therefore only the surface-to-surface performance comparison 
was made. 

The dispersion on target improved for the tubular 
projectile from 700 to 2000 meters. The dispersion varied 
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Figure 1. 
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20-rnm Ball Obturated Tubular 
components , and conventional 
(Photograph courtesy of NWC, 



Projectile (BOT) , 
20-ram projectile. 
China Lake) . 




Figure 2a. Standard 20-mra projectile in flight with 
detached bow shock. (Photograph courtesy 
of NWC, China Lake) . 
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Figure 2b. BOT with passage blocked resulting in 
detached bow shock. (Photograph 
courtesy of NWC, China Lake) . 




[ Figure 2c. 



BOT with no blockage in passage; 
oblique shocks emanating from the 
projectile. (Photograph courtesy 
Lake) . 



only attached 
lip of the 
of NIVC, China 
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on the average by only 0.2 mil over the entire range spectrum 
between the tubular and conventional HEI projectile. However 
at 2000 meters, the tubular was better by 0.4 mil. The time 
of flight to target was another impressive finding in Ref. 5. 
The tiobular projectile had a 30 percent shorter flight time 
at 1000 meters and a 40 percent shorter flight time at 2000 
meters than the conventional HEI projectile. 

Results of firings of 25-mm and 30-ram tubular pro- 
jectiles from the Oerlikin KBA and GAU-8/A respectively were 
reported in part in Ref. 3. Figures comparing Line-of-Sight 
Penetration Capability vs. Range for a conventional spinner 
configuration based on AR-2 shape, a f inner configuration 
based on FlO shape, and a tubular projectile were shown. The 
tubular projectile was shown to be inferior in both instances 
The report does not indicate the degree of compatibility the 
tubular configuration has with the guns used. 

4 . Review Of Internal Shock Wave Considerations 

From photographs of tubular projectiles in flight [4] 
and during wind tunnel testing [7], various shock patterns 
have been observed at the inlet region and in the wake. The 
photographs of Figs. 2a-c show the BOT and the conventional 
projectiles in flight. The sabots used in Refs. 4 and 7 were 
of the pusher type similar to those illustrated in Fig. 3. 

The detached bow shock seen in the photograph of the standard 
projectile (Fig. 2a) is similar to that of the tubular pro- 
jectile shown in Fig. 2b with passage closed (and in Ref. 4 



22 





23 



Figure 3. Solid (a) and Split (b) disk, pusher-type, discarding sabots/obturators. 



with the sabot attached) . With no blockage of the internal 
passage (Fig. 2c) there was no detached bow shock, with only 
attached oblique shocks emanating from the lip of the pro- 
jectile. This latter case allows the possibility of four 
other internal flow/shock configurations. 

(1) A normal shock standing at the entrance to the pro- 
jectile; 

(2) A normal shock standing at some intermediate position 
in the channel; 

(3) A system of oblique shocks present in the channel; 

(4) The channel is devoid of shocks and the flow is shock- 
free through the channel. 

In each of the cases described, specific conditions 
must exist; 

a. The Detached Bow Shock 

When the projectile passage is closed, such as 
in the case of the attached sabot, the projectile acts as a 
blunt-nosed body with the resultant detached normal shock (10] 
(see Figs. 2a and b) . From oblique shock theory, there is a 
maximum turning angle which the flow (supersonic) can nego- 
tiate through an attached oblique shock. When this turning 
angle is exceeded, the shock becomes normal and detaches. 

There have been many studies as to the strength and location 
of this detached bow shock, but here the internal flow is of 
more interest and, therefore, those references are not in- 
cluded in this discussion. 
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b. Normal Shock At The Entrance 

From one-dimensional frictionless flow theory, 
for a given initial Mach nvunber, , the pressure ratio 
across a normal shock is given by; 



P 

2 



P 



1 



(1 + ) 
1 



(1 + YM^) 
2 



( 1 ) 



for a perfect gas with constant specific heats. P is the 

2 

pressure downstream of the normal shock, P^ the pressure up- 
stream of the shock, M the downstream Mach number and y the 

2 

ratio of specific heats, Cp/C^. This pressure ratio must 
exist in order for a normal shock to be formed. Kantrowitz 
and Donaldson [11] conducted experiments in which they con- 
cluded that the contraction-ratio in the entrance of a super- 
sonic diffuser was important in reducing the kinetic energy 
losses. The studies done were on supersonic diffusers with 
varying areas. The formation of a normal shock at the en- 
trance is determined by the throat area and initial Mach 
number such that the Mach number at the throat of the duct 
is exactly 1.0. The maximum value of the ratio of the capture 
(inlet) area, A^, to throat area, A^ , required to allow the 
bow shock to be swallowed is given in Ref. 9 as: 

Y+1 

(2) 



_£ = l_ 

At M 



[ 



Y+1 



2(y-1) 



(1 + V ^2^J 
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where may be obtained from tabulated values for normal 
shock waves for a given free stream Mach number, M or may 
be calculated from 



2 

c. Normal Shock At An Intermediate Position 

The formation and position of a normal shock at 
an intermediate position in a constant area duct has been 
theorized in the presence and absence of viscous effects. 

In the absence of viscous effects, the noarmal shock tends 
to grow from the interaction of oblique shocks and will be 
discussed under "d" below. 

In the presence of friction, the formation of 
the boundary layer along the duct walls reduces the effective 
area of the duct [12, 13, 14]. This reduction has the effect 
of decelerating a supersonic flow as in a diffuser. In a 
constant area duct with friction, the location of a normal 
shock in the duct may be determined by Fanno Line Flow Theory. 
The length to diameter ratio, friction factor and initial 
Mach number determine the range of back pressures for which 
the duct will be "choked", or a normal shock set up [10]. 

For a given supersonic entrance Mach number, there 
exists a maximum value of fL/D; f, being the Reynolds number- 
dependent friction factor, L, the length of the duct and D, 
the diameter of the duct. This maximum value, is the value 



(Y-1) 



+ M" 



-L 



(Y-1) 



M"^ - 1 



(3) 
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for which the duct will remain unchoked. If the actual value 
of fL/D is less than the maximum, the flow will remain un- 
choked if the back pressure is sufficiently low at the duct 
exit. 

In the case of the tubular projectile, assuming 
a constant back pressure, L/D of (typically) 8.0 and Re in 
the range from approximately 4x10® to 6x10®, from Fig. 20.1 
of Ref. 15, f is between 0.013 and 0.0125. Using an average 
value of f of 0.01275, 

fL/D = 0.102 

From Ref. 10, Appendix E, this corresponds to M=1.41. 

Iterating to obtain better values of f from a revised average 
Re would increase the Mach number at which the projectile 
would choke. After one iteration the Mach number increases 
to M = 1.44. As L/D decreases, Re increases and f decreases 
resulting in a lower value of the Mach mamber at which the 
flow chokes. 

d. System Of Oblique Shocks In The Duct 

The formation of a system of oblique shocks in a 
duct or channel is the most difficult to explain and analyze. 
Ferri [16, 17] has developed an analytical solution for the 
shock structure utilizing the Method of Characteristics. 

The analysis becomes difficult if not impossible at the 
centerline of the duct because, according to Ferri, the oblique 
shock curves toward the centerline and becomes normal. The 
extent of the normal shock is dependent upon geometry and flow 
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conditions . The larger the wedge angle for a given Mach 
number, the stronger and larger the normal shock at the 
centerline. The constant area duct implies a zero internal 
deflection angle. In this instance, Ferri ' s analysis, which 
ignores viscous effects, implies that no normal shock is 
formed (Fig. 13 of Ref. [16]) . With no wedge angle, there 
is no oblique shock formed. 

Bruns void and Kalivretenos [9] have modified a 
Method of Characteristics computer code [18] developed for 
supersonic inlets to account for the absence of a center-body 
on the tubular projectile. This modified code was then 
coupled with a subroutine to calculate the internal wave drag 
coefficient based on and used to obtain an optimum internal 
geometry for the tubular projectile which would give minimxam 
internal wave drag. The computer code assumes inviscid, 
supersonic flow throughout. Reference 9 states that some 
inaccuracies are allowed in the computation due to the re- 
laxation of the convergence of the characteristics equations 
at the centerline, but that the inaccuracies do not have time 
to propagate to the walls of the projectile before the exit 
plane. Therefore, the inaccuracies do not affect the wall 
pressure distribution. 

Rhethorst, et.al. [19] has postulated another 
method for determining the shock structure and the resultant 
internal wave drag for a tubular projectile, using Prandtl- 
Meyer wave theory. But again the initial oblique shock 

2 8 



formation is dependent upon the wedge angle at the entrance 
and therefore does not apply to the constant area duct case 
where no internal wedge angle exists. Ferri [16, 17], 
Brunsvold and Kalivretenos [9] and Rhethorst, et.al. [19] 
have neglected boundary layer/shock wave interaction. 

Fejer, et.al., [14] and Waltrup and Billig [15] 
have conducted experiments utilizing constant area ducts. 

Fejer, et.al., performed experiments on constant area rec- 
tangular ducts varying in lengths, whereas Waltrup and Billig 
conducted their experiments on cylindrical ducts of varying 
lengths and diameters. Other than the obvious geometry dif- 
ferences, the procedures were very similar. The advantage 
of the rectangular duct was the ability to visualize the shock 
patterns through a viewing port, whereas the curvature of the 
walls of the cylindrical duct precluded accurate visualization. 

In both experiments, oblique shock systems were 
formed in the ducts. The location of the shock system in the 
duct and the actual configuration of the shock system was con- 
trolled by the pressure ratio. The pressure ratios were con- 
trolled in different ways, but the effects were the same. As 
back pressure increased, a shock system formed at the exit of 
the duct. With further increases in back pressure the system 
of shocks moved upstream and out of the duct and a detached 
shock appeared. Shapiro [12] indicates that internal 
oblique shocks may be the remnants of normal shocks "with 
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bifurcated ends as a result of boundary layer separation." 

He states that the normal portion grows shorter as the 
boundary layer gets thicker until the normal portion dis- 
appears completely leaving only the bifurcated portion. He 
illustrates this in Fig. 28.27(c) of Ref. 12. This explains 
the movement of the oblique system toward the entrance of the 
duct in the experiments of Refs. 14 and 15. As back pressure 
increases, the flow downstream of an internal shock has more 
of a tendency to separate. Fejer also states that the 
presence of a turbulent boundary layer had a stabilizing 
effect on the shock system in the constant area rectangular 
duct. In the experiments involving the cylindrical duct, the 
shock structure again moved from exit toward entrance with 
back pressure increases. They also report a breakdown in flow 
downstream of the shock system as indicated by relatively 
constant pressure movements in this region. 

e. No Shocks Present In The Duct 

It can be seen from previous discussion that 
given an initial Mach number and inlet stagnation pressure, 
a sufficiently high back pressure will cause a normal shock 
or system of oblique shocks to be formed. If the cross- 
sectional area of the duct is reduced in the direction of 
flow, as in a supersonic diffuser, the incident flow may 
cause an oblique compression wave to be formed. In the 
absence of the requisite back pressure or area reduction. 
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no explanations for internal shock formation have been fo\ind 
other than those due to viscous effects noted in "d" above, 
f . Summary 

In a constant area duct with supersonic inlet 
flow, there is a pressure rise caused by a reduction in the 
effective flow area by boundary layer formation due to 
friction. Depending upon the length-to-diameter ratio (L/D) 
of the duct and the flow parameters, the increase in pressure 
ratio (or reduction in Mach nximber) required to cause forma- 
tion of a system of oblique shocks is less than that required 
to set up a normal shock at the duct entrance. This difference 
is explained by the boundary layer formation and/or interaction 
with the shock waves. 

In the case of the tubular projectile, the back 
pressure at the exit in flight is affected by the complex wake 
flow. There is no pressure difference in the surrounding 
medium at large distances from the projectile. The formation 
of shock waves inside the projectile with a constant cross 
section is due to viscous effects and depends upon Mach number 
and the length-to-diameter ratio. 

5 . Laboratory /Wind Ttinnel Testing 

From all indications, the Weapons System Concept Team 
(WSCT) of ARRADCOM has done extensive wind tunnel 
testing on various tubular projectile configurations. The 
work of A. Plateau of the Weapons System Concept Office is 
referenced often with regard to results obtained, however no 
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publications from that office have been found. (The work 
referenced in other publications is usually of the "private 
communication" type.) Appendix D of Ref. 8 reports results 
of a test conducted by the WSCT at the ARRADCOM facility. 

The facility utilizes a variable Mach number tunnel which 
may be used to simulate the deceleration of the projectile 
in flight. 

The results of Ref. 8 agree closely to those of 
Ref. 7, in that the drag coefficient rises almost instan- 
taneously at the Mach number at which the flow becomes 
choked. 

The Mach number at which this choking occurs is de- 
pendent upon geometry. For the case of inviscid flow, 
Brunsvold and Kalivretenos [9] give the following relation 
for the area ratio corresponding to the free stream Mach 
number at which the shock is "regurgitated" and the flow be- 
comes choked: 





Y+1 

_,2(y-1) 




(4) 



C. PURPOSE OF STUDY 

The results of tests and experiments previously mentioned 
show for the most part definite advantages to the tubular 
projectile over the conventional projectile primarily in the 
area of lower drag (approximately 1/2 to 1/3 of conventional [8]) 
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which allows a flatter trajectory and shorter time of flight 
for a given range [5j . With the emphasis on Anti-Missile 
Defense both from a Surface-to-Air and Air-to-Air standpoint, 
these features are definitely worth further study. 

The problem seen from the air platform is that of the 
discarding sabot or pusher. The solid [4, 5, 7] or split [8] 
disk is effective and reliable, however the possibility of 
ingestion in the aircraft engines prohibits their use. 
Rhethorst et.al. in Ref. 7 illustrates a number of sabot/ 
obturator designs for use with the tubular projectile, but 
these are also of the discarding type. They also investigated 
the possibility of using a consumable sabot [19] which would 
be burned up as the projectile was transiting the barrel and 
would be completely consumed as it exited the barrel. The 
most recent development has been the EOT which was designed 
at NWC China Lake, California [6]. The ball obturator has 
been bored with a hole the same diameter as the hole through 
the projectile. When loaded, the ball is supported inside 
the projectile such that the axis through the hole in the ball 
makes an angle of approximately 90° with the axis of the hole 
through the projectile (see Fig. 4) . 

Gas pressure from the burning propellant holds the ball 
obturator fixed with respect to the projectile as the rifling 
in the barrel imparts spin to the projectile (as well as the 
ball) . Upon exiting the barrel, the propellant gas pressure 
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Figure 4. BOT with obturator oriented to block the passage (loading position) . 



is released. At that time a complex force distribution acts 
on the ball to cause it to change position inside the pro- 
jectile in such a way as to align the holes. Inertial forces 
then dominate to maintain this alignment so that the pro- 
jectile remains fully tvibular. There are no separating parts 
and the opening process is automatic. 

The purpose of this research is to analytically predict 
the motion of the ball obturator as a function of time. In 
support of this purpose, an experimental apparatus was de- 
signed and built to simulate the motion of the obturator in 
order to gain physical evidence of the accuracy of the model. 
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II. ANALYSIS 



A. DEVELOPMENT AND SOLUTION OF THE EQUATIONS OF MOTION 

The ball obturator may be thought of as a rigid body of 
rotation with a system of coordinate axes fixed to the ball 
having its origin at the mass center of the ball (Fig. 5) . 

These axes will be designated x, y, and z, where z is the 
axis through the hole in the ball. 

The motion of the ball must be described relative to the 
projectile, and for this purpose an inertial frame of refer- 
ence is defined with axes designated X, Y and Z with its 
origin also at the mass center of the ball. For this analysis 
it is assumed that the projectile has only a rotation about 
the Z-axis and no translation or other rotation. Because only 
the motion of the ball with respect to the projectile is de- 
sired at this point, the relative velocities are the quantities 
desired. The ball has only rotational components and thus the 
values of angular velocity with respect to the X and Y axes 
are also relative to the projectile, whereas the value with 
respect to the Z-axis must account for the projectile spin. 

By doing this, the motion may be described by Euler's 
Modified Equations of Motion for a rigid body of revolution 
about a fixed point. These are given in Thomson [20] as: 
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Figure 5, BOT showing orientation of coordinate axes. 



(5a) 



A(I) + (C-A) 0) w = G 

1 2 3 1 

Ad) + (A-C)o) w = G (5b) 

2 13 2 

Cd) = G (5c) 

3 3 



where 

A = Mass moments of inertia about principal 

axes perpendicular to the z-axis. 

C = Mass moment of inertia about the prin- 

cipal axis through the hole (z-axis) . 

(I) ,d) ,d) = Angular accelerations about the body 

^ ^ ^ fixed axes x, y, z respectively. 



(i) / (l) r 00 

1 2 3 



Angular velocities about the body fixed 
axes . 



G ,G ,G 

1 2 3 



Moments acting about the body fixed axes. 



Since a body of revolution has only two distinct prin- 
cipal axes, the mass moments of inertia about the x and y 
axes are equal. The moments of inertia of the ball were 
determined by first finding the inertia tensor. From 
Refs. 21 and 22 the inertia tensor is found by starting with 
the formulation for a hollow cylinder of inside radius r and 
length 2h (Fig. 6) in which the inertia tensor in terms of 
dyadics is 

I = (U + ee)mr^/2 + (U-ee)mh^/3 (6) 



In Eq. (6) m is the mass of the cylinder (ball) , U is the 
unit dyadic and ^ is the referenced dyad. Using Fig. 6 as a 
guide, the incremental value of the inertia tensor with respect 
to the center of mass, dl^^ is given by: 
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Figure 6 



( h = /R^-p2 ) 




Dimensional relationships for calculating the 
inertia tensor and moments of inertia. 



.39 



dip = (U + ee)p^/2 dm + (U - ee)h^/3 dm 



(7) 



where p is the radial distance to the incremental mass, dm. 
The incremental mass is given by 



dm = 4irp(R^~p^) ydP 



( 8 ) 



From Eq. (7) , since the origin and center of mass coincide. 






dl. 



(9a) 



therefore , 



I = 
o 



rR u 3 

/ (U+ee) 2p ^TThydp + / (U-ee) j— 4iTpydp 
Jr Jr 



(9b) 



which, after simplification leaves 



rR h 

I = 27ry(U+ee) p^R^-p^) dp + 



(9c) 



■R 



3/2 



•|TTy(U-e£)| (R^-p^) pdp 



Solving Eq. (9c) gives 



3/2 



I = iry{2(U+ee) [ir^+i-RMR^-r^) ] + 

o — b 15 



4 rl 2 2 

§(U-ee) [j(R2-r^) ]> 



(9d) 
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Eq. (9d) is the expression for the second order inertia 



tensor of the ball obturator with outside radius R, hole 
radius r and made of a material of density y. Diadics have 
been used to describe this tensor. In order to obtain the 
moments of inertia about any axis, the diadic describing the 
tensor need only be premultiplied and postmultiplied by scalar 
multiplication by the unit vector along the axis about which 
the moment of inertia is described. It may be seen in Fig. 6 
that the e-direction vector is along the z-axis. The 
e-direction may be chosen arbitrarily, however this choice 
greatly simplifies the formulation of the integral. Thus 
the moment of inertia about the z-axis, C, is given by 



Simplifying and using the subscript s to denote particular 
values for the ball obturator: 



For the axes perpenducular to e through the mass center, the 
Moment of Inertia is given by 



1 2 

C = e • [ (ii+j j+ee) +ee) ] * e { [2 t7y ( jr^+j^R^ (R^-r^) 



3/2 



4 3^ J/ ^ 

+ -(ii+j j+ee-ee) [iry-g- (R^-r^ ) ]} 



(10a) 



3/2 




(10b) 




(11a) 
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Performing steps as in Eqs. (10) gives 



3/2 

A = j • (ii+j j+ee+ee) • j [2ity (^r (R^-r^ ) ] 



5/2 

+ j • (ii+j j+ee-ee) * j [Try-ijCR^-r^ ) ] 



(11b) 



After premultiplying and postmultiplying, 

12 3/24 5/2 

A= iTY{2[p:s^ + ^(llc) 



In order to simplify Eqs. (5) , define a non-dimensional 
X such that 

X = ^ (12) 



Thus Eqs. (5) becomes 




(I) + Xcow = G/A 

1 2 3 1 


(5a) 


u) +Xcow = G/A 
2 13 2 


(5b) 


d) = G /C 

3 3 


(5c) 



The position of the ball at any time may be described by 
three angles, 9, (J) and ip. These are known as the Euler angles 
and their relationships to the coordinate axes already de- 
scribed may be seen in Fig. 7. Reference 20 defines 

and 0 ) in terms of the Euler angles as 
3 
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Figure 7. Relationship of Euler Angles, 0, i|; and (j) to the coordinate 
axes. (c is the nodal axis in the X-Y plane). 
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(13a) 



0 )^ = ^sin6sin(j) + 6cos(j) 

= iJ)sin6cos<j) - 0sinc{) 

• 

= (}) + ij;cos6 

By differentiating Eqs. (13) with respect to time, w , 

1 

and w may be obtained as 

3 

03^ = i|)sin0sin({) + ( 0COS 0sin({) + (f)sin0cos(j)) + Bcos<{) 

• • 

-0(j)sin(i) 

• •• • • • •• 

03^ = ii/sin6cos(() + (0 cos0 cos<j> - 4>sin0sin<|) ) - 0sin4> 

- 0^COS(}) 

• • 

J) = $ + i^cos0 - i|)0sin0 

3 

By substituting Eqs. (13) and Eqs. (14) into Eqs. 
the following expressions for G /A, G /A and G /C are 

12 3 

to be 

• • 

i|}sin0sin({) + i|3 (0 cos 0s in <{) + 4>sin0cos(j)) + 0cosc{» - 
0^sin({) + X (^sin0cos(j) - 0sin(f») (ii3cos0 + ^) = G^/A 

ipsin0cos(J) + 1^3 ( 0cos0cos4) - ^sin0sinc{)) - 0sin(j) - 
0^cos(J> - X (ii3sin0sin(i) + 0cos4>) (4»cos0 + ^) = G^/A 

$cos0 - Tj50sin0 + $ = G /C 

3 



(13b) 

(13c) 

w , 

2 

(14a) 

(14b) 

(14c) 

(5) , 
found 

(15a) 

(15b) 

(15c) 



Transforming this system of differential equations into the 
matrix form [A] {x} = [b] leaves 



m 











— 








( y 




sin9sin(j) 


0 


cos4> 




$ 




sin9cos(j) 


0 


-sincf) 


< 


$ 


> = 


cos 9 


1 


0 




0 












V ) 


L 



G^/A - ^ (6cos0sin(J) + ^sin9cos(f)) 

• • • 

G /A - ij; (6cos9cos(J) - (j)sin9sin4)) 
2 

G^/C + ^0sin9 



+ 0^sin4) - X {^sin9cos(}) - 9sin(I») (^cos9 + 4>) 
+ 0^cos(j) + X {ij;sin9sin4) + 6cos$) (iijcos9 + 



(16) 



By premultiplying both sides of Eq. (16) by the inverse 
of the coefficient matrix, we are left with 







sin4> 


COS(f) 


0 




sin9 


sin9 


u 






sin4>cos9 

sin9 


-cos9cos(|) 

sin9 


1 


s 




COS(|) 


-sin<}> 


0 


/ 









G^/A - iji ( Scos 9sin4> + (psin9cos<j)) 

• • • 

G /A - if) (9 cos 6 cos (j) - (j)sin9sin(j) ) 
2 

G /C + ^9sin<|) 

3 



+ &4>sin<}) - X (if)sin9coS({) - 9sin(})) (ificos9 + ^) 
+ 9^cos4> + X (i|)sin9sin4) + 0cos4>) (iJ)COs9 + ^) 



(17) 



In order to solve this system of non-linear second order, 
differential equations for 9, (j) and \p, the following must be 
known ; 

• • 

Initial Values (denoted by subscript o) of 9, <p, ijj, 9, <p 
and ij; and the quantities. A, C, G , G and G . 

12 3 
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Due to symmetry, 0° £ 6 £ 90°, and since we are interested 
in the performance of the ball away from 0°, 9^ > 0. Since 

the ball is attached to the projectile initially (no slippage) , 

• • • 

0Q = 4>q = 0. equal to the spin rate of the projectile 

as it exits the barrel. A and C are physical constants de- 
pendent upon geometry and material properties [Eqs. (10) and 
(11) ] and therefore X is a constant. <t>^, and are arbitrary 
due to symmetry and the values of G , G and G must be de- 

1 2 3 

termined from knowledge of forces acting on the ball. 

The solution to the problem of the motion of an axi- 
symmetric rigid body about a fixed point with no moments 
applied is well known [20,21,22,23]. Under these conditions 
there is no nutation and the precession rate ^ is constant. 

Thus 

6 = 0 



and 



ip = 



C(j> 

(A-C) COS0 



(18) 



The computer solution method was tested by substituting the 

projectile spin velocity for ip in Eq. (18) and calculating 

• o 

(}) for the specified value of 0 < 6 q < 90 . By inserting 

these values of 0^, i)) and ^ into the computer model, and 
solving for 0 as a function of time, 0 was indeed found to 

be zero. In addition, all other criteria for the test case, 

• • 

i.e., constant angular momentum and constant \p and <j) were met. 
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The problem remaining is one of modeling the forces acting 
on the ball in order to determine G , G and G . 

12 3 

B. MODELING OF THE APPLIED TORQUES 

1. Moments Due To Viscous Effects 

Without the influence of applied forces or moments, 
the angular momentum remains constant or, from first prin- 
cipals : 

t 

M = h = 0 (19) 

In order for the ball to nutate to the position in 
which the holes are aligned, a moment must be applied. This 
resultant torque is due to both viscous and sliding frictional 
forces . 

The first attempt was to model the torque on the ball 
as that due to laminar Couette flow of a viscous fluid between 
two flat plates where the shear stress is given by 

T = (20) 

The solution of the Navier-Stokes equations for this 
flow is well known and given by 



dP 

dx 




( 21 ) 



where dP/dx is the pressure gradient in the direction of 
flow and, for our calculations, is assumed to be zero. The 
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governing partial differential equation is then 




0 



( 22 ) 



By considering the relative velocity, one plate is held 
fixed while the other has velocity of magnitude w^r. coj. 
is the relative angular velocity about the instantaneous 
axis of rotation and r is the perpendicular distance from 
the instantaneous axis to the point at which the velocity 
is desired. 



where y=0 is the surface of the ball and h is the gap width 
between ball and projectile. Integrating Eq. (22) twice and 
applying the boundary conditions Eq. (23) gives 



Thus from Eqs . (20) and (24), the viscous shear 

stress at the surface of the ball, or at y=0 becomes 



The boundary conditions are therefore given as 



u = -ojj.r at y = 0 



(23) 



u = 0 



at y = h 



u(y) = -Wj.r (1-y/h) 



(24) 



0) r 



T 




(25) 
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The force due to viscous shear on the ball is given 
by the integral over the surface area of the ball of the 
shear stress Eq. (25) , or 



f'v = I "o“- 



(26) 



The differential area dA is given by 



dA = 2 ttR_ (R_sin(xda) 



(27) 



which is shown in Fig. s* 

Substituting Eqs. (25) and (27) into Eq. (26) gives 



Fv = 



2ru^Rs"Wr 



IT 

rsinada 



(28) 



The value of r changes from 0 at the axis of rotation 
to Rg at a=90°. Thus 



r = RgSina 

Substituting this expression into Eq. (28) 



(29) 



F 



V 



2ttu.R ^oj 
f s r 



h 




IT 

sin^adct 



(30) 



o 

The torque imparted to the ball is the force given by 
Eq. (30) multiplied by the moment arm at which it acts. 
Again, the moment arm is a function of Rg and a and is given 
by RgSina. Substituting this into Eq. (30) gives an ex- 
pression for the torque acting on the ball obturator due to 
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Figure 8. Differential area, dA, orientation with respect 
to the relative angular velocity vector, 
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the velocity of the ball relative to the projectile in a 
viscous medium, or 



27ry R r 

Q = / sin^ada (31) 

o 

Carrying out this integration gives 




TTRs'^cOryf 

h 



(32) 



For a given geometry and constant fluid properties, 
co^ is the only variable in Eq. (32), therefore that portion 
of the expression which is multiplied by will be called 
the shear constant, S, such that 

Q = SWj. (33) 



where 



S 





(34) 



From the rigid body motion analysis, the relative 
angular velocity is easily determined. It is 

u = [(u -w )^ + (u -CO )^ + (CO -CO )^] (35) 

r 1 Pi 2 P 2 3 P3 

where co_ , co_ and co„ are the components of the projectile 
Pi P2 P3 

angular velocity about the axes fixed in the ball. These 
relationships are given in Thomson [20 ] 
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(36a) 



w = 03^sin(|>sin9 
Pi P 

Wp^ = WpCOS(j)sin0 (36b) 

01 = 0 ) COS0 (36c) 

p3 p 

Equations (36) are simplified by the fact that oip 
has only a component in the Z-direction of the inertial 
frame of reference. 

The moments G , G and G are obtained by applying 

12 3 

Eq. (33) and resolving the moment vector into components 
about the x, y and z-axes. The motion of the ball may then 
be described fully by Eq. (17). 

Figure 9 shows solutions for 0 as a function of time, 
using the Couette flow model with the following inputs: 

R = 7.62 nun (0.3 in.) 

o 

rg = 4.7625 mm (0.1875 in.) 

h = 0.0254 mm (0.001 in.) 

p = 7800 Kg/m^ (0.282 Ibm/in^) (steel) 

= 1.92 10“® Pa-sec (4xl0“ ^Ibf-sec/ftM 
Wp = 12042.77 rad/sec (1. 15 x1o5rpm) 

It appears that viscous effects alone are not suffi- 
cient to cause the ball to nutate toward 0=0 in a rapid 
fashion. To check the model again, progressively larger 
orders of magnitude of S were input into the model. The 
effect, as expected, was to cause 0 to approach zero more 
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rapidly Csee Fig. 10). From firings conducted at NWC, China 
Lake, on an’ indoor range, it is known that 9 reduces to a 
value which allows light from behind the projectile to be 
seen through the projectile by a camera aiming at the pro- 
jectile image in a mirror down range. The length of the 
range is such that, for this to occur, 0 must reach the value 
3g prior to impact on target, or in about 0.1 seconds. (see 
Fig. 11 for 3g) . The viscous model described above does not 
produce the moments required to cause sufficient reduction of 
9 within this time of flight. 

The viscous torques on the obturator were calculated 
above assxming a solid sphere. The actual obturator has an 
outside surface area less than that of a solid sphere of the 
same radius. The surface under the influence of viscous 
shearing stress is a very complicated function dependent upon 
the angle 0 as shown by the cross-hatched portion of Fig. H. 
Calculations of the resulting moment are further complicated 
by the relative effect of the missing areas (due to the 
presence of holes) and their positions relative to the 
instantaneous axis of rotation. The time-dependent angles 
between the instantaneous axis of rotation and the z and 
Z-axes must be known in order to determine the limits of 
integration for correct moment arm consideration. 

Because the moment due to viscous effects is dependent 
upon the surface area involved and because the surface area 
in the first calculation (assxaming a solid sphere) is larger. 
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it is to be expected that by taking the holes into account 
the true moment due to viscous effects will be reduced. The 
effects of these area variations will be "blurred" by the 
high-speed relative motion and the additional complications 
necessary to account for them are not warranted in light of 
the goals of this analysis. Therefore for further develop- 
ment of viscous effects, the surface area of the obturator 
will be considered that of a solid sphere. 

Nakabayashi [24] has conducted experiments to de- 
termine the viscous torque on rotating concentric spheres 
with various fluids and gap widths. He gives the general 
expression for the frictional moment on a sphere due to flow 
in the gap as 



M 



V 







(37) 



where p is the density of the fluid in the gap. For laminar 
flow 



C 

m 




(l+£) 

e 



Re"^ 



(38) 



and for the turbulent flow 

Re “0.25 

„ _ 0.053 l+(7/4)e 

m e 1-{3/2)£ 



(39) 



where e = h/Rg. 

The Reynolds number is defined as 

R 

Re = — — ^ (40) 

V 
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« 

m 







in which v is the kinematic viscosity of the fluid in the 
gap. For the narrow gap width assumed for the BOT, the 
transition Reynolds number may be approximated by 

Re^ = 70e"l*^ (41) 

(from Fig. 2a of Ref. 24.) in the case of the BOT, the 
relative angular velocities are insufficient to assxame the 
existence of turbulent flow in the gap. Thus the frictional 
moment due to viscous flow in the gap is given by 



Mv = - 




(l+£) 

e 








(42) 



For small values of e (as is the case here) , Eqs. (42) and 
(33) are identical. 

2 . Moments Due To Sliding Friction 

It has been noted that the extreme propellant gas 
pressures inside the barrel are assumed sufficient to hold 
the obturator fixed to the projectile during launch. How- 
ever upon exiting the barrel/ the pressure distribution on 
the projectile is due to the aerodynamic characteristics of 
the projectile. Assuming the obturator is still blocking 
the projectile duct, a detached bow shock is set up (see 

Fig. 2b) . The pressure, P , downstream of this shock is given 

2 

by one-dimensional shock theory from Eq. (1) . It is assxamed 
that is sufficiently large to cause the obturator to con- 
tact the rear of the projectile socket housing the obturator. 
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The forces involved depend upon the projected area 
of the hole through the projectile and the contact point (s) 
of the obturator and projectile. The force acting on the 
obturator is assumed to be given by 

F = irr 2 AP [see Eqs. (1) and (3)] (43) 

P P 

where r is the radius of the hole through the projectile 
and 

AP = P - P (44) 

1 2 

From sliding friction theory, the contact area is 
not a consideration; only the magnitude of the force normal 
to the contact surface and the coefficient of friction of 
the materials in contact are important. In the case of the 
obturator, the pressure force, F acts along the Z-axis. 

The reaction force may be broken into a component tangent 
to the obturator and a normal component acting through the 
center of the obturator. The normal component, F^, multi- 
plied by the coefficient of friction Pg, gives the sliding 
friction force acting to oppose the motion of the obturator; 
Fg = ^n^s 

where 

= FpCosSg (45b) 

so that 

F^ = F^cos3^U^ 
s p s s 
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This assumes the holes through the obturator and projectile 
have the same diameter and that the point or points of contact 
are on the locus of points described by the intersection of 
the hole through the projectile and the spherical socket 
housing the obturator in the projectile. This locus of points 
is a circle of radius rp and each point is a distance Rg from 
the center of the obturator. Gravity forces are neglected 
since they are relatively small and act both on the obturator 
and the projectile. 

For a "perfect" fit between ball and projectile, 
the sliding friction force is distributed over an infinite 
number of points and, in the limit, the force at each point 
vanishes. This situation has no significance in reality and, 
therefore, engineering judgement dictates making some assump- 
tions with regard to the point (s) of application of the 
friction force Fg. 

For the general case, a point of contact P may be 
described relative to the inertial reference frame as being 
at an angle 3g from the inertial Z-axis, a distance Rg from 
the center of the obturator and at an angle y from the 
X-Z plane measured toward the positive Y-axis as shown in 
Fig. 12 . For this approximation, assume two contact points 

in the X-Z plane, and P^. P^ and P^ are 180° apart on 

the circle of points so that y=0 and 6=Sg. 
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From the vector analysis, the projection of one vector 
onto another is given by the scalar or dot product. The 
vectors and are given by; 



P = Re (sinBi + cosBk) 

1 s — — 

P = R (-sinBi + cosBk) 

2 s — — 



(47) 



Their projections onto the relative angular velocity 
vector ( 0 ^/ with components Uy and are 



P *a)_ = R_ (co^sinB + o)„cosB) 

1 ^ X 2 

(48) 



P *03^ = R^ (-w^sinB + o),cosB) 

2 i S X Z 



The scalar dot product is also defined as 



P = P w^cosri (49) 

1 ^ 1 ^ 

where n is the angle between the vectors P^ and co^^. Therefore 
from Eqs. (48) and (49) 

R (o) sinB + 0 ) cosB) = P o) cosn 
s X z 1 r 

(50) 

R^(-o)„sinB + o),cosB) = P u^cosfi 

S X Z 2 ^ 



The magnitude of the vectors P and P is R , so solving 

1 2 S 

for n and Q gives 
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03 03^ 

Ti = cos"*^ — sing + — cos6 

03 03 

r r 



( 51 ) 



Q = cos""^ - — sing + — cosg 



03, 









CO, 



The moment arms and from Fig. 12 are therefore given by 



03 03 

. - 1 / X^,. . Z. 



r = Rc.sin [cos” ^ (-Rising + «^os3) ] 



(0, 



03. 



03 03 ^ 

r = R [1 - (-^in6 + -^osB)^] 

1 S 03 03 



(52) 



03 03 

X . ^ Z 



r = R sin [cos"^ (- ^=>xnp -r 



inB + e^cosB) ] 



r = R [1 - (^^osB - ^^inB)^] 

2 S (1)^ 0)^ 



(53) 



The moment acting on the obturator due to the sliding 
friction due to aerodynamic forces are given by 



->• 

03 



M = (r +r ) ^ cosBy- — 
Sa 12 2 s 



(54) 



Substituting from Eqs. (52) and (53) into Eq. (54) leaves 



M, 



_ /2 



sa 



R FpCosBy_[l - (— )^sin^B - (— )^cos^B]^ ^ 

S h' S U\ (jy 



“2 



w. 



(55) 



S2 



It should be noted that if there were n contact points that 
the normal force at each point would be given by F^/n. 

Equation C54) indicates the case of two points of contact 
discussed here. 

As yet, there has been no attempt in the laboratory 
to produce the aerodynamic forces present on the obturator 
for an actual projectile in flight. The sliding friction 
forces present on the obturator in the experiments (described 
in section III) are due to gravity alone. The obturator is 
supported by the projectile simulator (see Fig. 24 )• Because 
the obturator is smaller than the housing in the projectile, 
a point contact at the bottom of the housing is assumed. 

This point has inertial coordinates of (“Rg, 0, 0). The 
vector from the origin to the point is given by 

->■ 

Pg = -Rsi (56) 

Using a development similar to that for obtaining the moment 
arm for the friction forces due to aerodynamic considerations, 
the forces on the model due to gravity may be determined. 
Again, the dot product is given by 

^G*“r ~ ~^s“x (57) 

which is the same as PqW^cosC • 

Therefore the angle (?) between the angular velocity 
vector Uj. and Pq is obtained from 




■m 




( 58 ) 



r 

The moment arm is found as in Eq, (55) from 

r = R sin [cos”^ (- ^) ] (59) 



or 



0) h 

= R3 [1 - (5^>n 



(60) 



which means the sliding friction moment due to the weight W is 

CO.. h CD. 



M3^ = y,WR, (1 - (j^) M 



(61) 



The torques acting on the obturator are the 
summation of the torques due to friction, both viscous and 
sliding (aerodynamic or gravitational) . The general form 
of the torques acting to oppose the motion of the obturator 
in terms of total moment and the unit vector 1^ along is 

G = “M^lj. (62) 

which, when broken into the various contributors becomes, 
for the projectile in flight, 

Gp = -(M^+Mga)lj. (63) 
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TABLE I 





COMMON 


VALUES 


So 


= 87.25° 




io 


= 0 




9o 


= 0 




<^o 


= 0 






= 0 




P 


= 7669.6 Kg/m^ 




Rs 


= 7.899 mm 






= 4.72 mm 




c 


= 3.125x10“'^ N-m-s^ 




A 


= 2.216x10”^ N-m-s^ 




X 


= 0.41 






= 1.917x10”^ Pa-s 






= 0.35 




h 


= 5.08xl0~® m 




S 


= 12.38x10“® N-m-s 




AERODYNA^lIC MODEL 


GRAVITATIONAL MODEL 




= 9.96 mm 


W = 0.07987 N 


"p 


= 76.2 mm 


E(6o,7r/2) = 1.0063 


m 

P 


= 0.141 Kg 


u = 200 sec * 

P 


I 

P 


= 8.58x10“® N-m-s^ 




M 

1 


= 3.0 




Y 


= 1.4 




P 

1 


= 0.099 MPa 




w_ 


= 12042 sec ' 





P 
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and for the experimental model. 



These may easily be transformed into the body-fixed coor- 
dinate system for solution of the equations of motion. 

For solving Eq. (17) , a numerical integration s\ab- 
routine [25] was used and computer plots of 9 vs. time were 
obtained. A copy of the calling program is included in 
Appendix A. Figures 13 and 14 show plotted output for the 
input parameters shown in Table 1. 



C. DEVELOPMENT OF AN APPROXIMATE SOLUTION 



• Observations of the nature of the exact solution (Figs. 

13 and 14) indicate that the response of the obturator orien- 
tation (0) to the applied torques (determined from Table 1 
parameters) is similar to that of a linear damped second-order 
system with a combined step and ramp input. That is, the 
relationship closely follows a relationship of the form 



9 



0O+D+ (w t-2?) 
CO n 
n 



^ exp(-co)i^t) 
( 1 -?^)^ 



{^in[to (1-^2) 't + (j) 1 
ti)„ n x\ 



u 

- D sin (1-?^ ) t + 4)p] } (65) 

in which the parameters are ; 

D = step height (negative) 

R = ramp slope (negative) 

= natural frequency 
C = damping ratio 
= phase angles 
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Figure 13. Plot of 0 vs. time for EOT utilizing aerodynamic force 
development. 
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Figure 14. Plot of 0 vs. time for EOT utilizing gravitational force 
development. 




a* 




These may be determined by means of a solution to the 
linearized equations of motion under the assximption of 
small applied torques [(M/Cw ^)<<1]. This has been accom- 
plished^ by a perturbation of the steady state solution given 
by Eq. (18) for the case of no applied torques. 

The first step in the approximation is to rewrite Eq. (18) 
to read 

COS0 + (X+1)^ = 0 (66) 

’^o o ^o 

• • 

where ({>^ and 0^ are constants. The dependent variables 

of the general form of the governing equations (Eq. 17) are 
written as sums of the unforced values and the perturbation 
values . Thus 

(67) 

$=$, ip=i)), 0 = 0, 9=0 

The general form for the applied torques on the obturator 
in the experimental model (gravity-driven sliding friction) 
is given by Eq. (64) in which the moment term is composed of 
a fluid (M^) and a sliding friction (Mgg) term. When both 
are expanded from Eqs. (42) and (61), Eq. (62) becomes 

^ w h to 

G = - {Sto^ + ygWRg [1-(^)^1 } ^ (68) 

r r 



2 

The theoretical development of the linear approximation 
is not complete. The analysis presented here therefore 
represents a preliminary report. 
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The relative velocity cj^ from Eqs . (13, 35 and 36) has com- 

ponents 

CO = (i-to )sin9sin(i) + 0cos<j) 
ir ^ p 

0),^ = sin9coS(J) - 9sin(j) 

^ (69) 

w ^ = (ij;-w„)cos9 + ^ 

3 r hr 

• • 

To the first approximation, ('Pq-^ the case of 

P P 

no applied torques) . With this approximation, and the 
application of the perturbed variables, Eq. (67), 

• 

oj = 9coS(l) 
ir - o 

cj = -9sin({) 

2r o 

• • • 

cj = i/;cos9,. + 4)_ + 4» (70) 

3r - u ^ 

The only term of zeroeth order in Eqs. (70) is (})q so 
that, to the first approximation, 

03^ = 0)^^ li^l (71) 

The modulus is taken in order to preserve the positive sense 
of co^. Under the same approximations as above, 

‘^x " ^^sin'i’^sin’i^Q (72) 

Substituting Eqs. (71 and 72) into Eq. (60) for the 
sliding friction moment arm, i^g/^s becomes 
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m 

% 



m 




( 73 ) 



r % 

^ = [1- (sine^sinil/Q) ^ J 

s 

For the purpose of the linear approximation, the mean value 
of this moment arm is required. Integrating to find the 
average moment arm gives 



- 

^ = i r ( 

Rs T J 



r /R ) dt 
g s 



(74) 



where T is the period for one revolution. From the approxi- 
mation, however. 



rp - 2ir _ 2TTt 
^o o 



(75) 



and because of symmetry, the average is valid over one-quarter 
period, therefore 



rt/4 



r 



4 

T 



2irt ^ 
(l-sin^6^sin^ — ^ ) dt 



(76) 



By letting u = 2TTt/T and du = (2iT/T)dt, Eq. (76) becomes 




_2 

IT 



7T/2 



/ 



o 



h 

(l-sin^6 sin^u) du 
o 



(77) 
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Letting E(0^) denote a Complete Elliptic Integral of the 
Second Kind with modular angle 9^, Eq. (77) may now be 
written as 



(78) 

Solutions to Eq. (78) are tabulated for various values of 
6o- 

The expression for the applied torque vector (Eq. 68) 
may now be written as 

G = -ISlil + UsWRg I 

With these approximations, the applied torque terms 
appearing in Eq. (17) are given by 



G G 

1 2 -M 

+ ^-cos(J) = jjj-(ii»-a)p) sin0-O 



%os4. - %in* = :^6 
^ “r a|4.„ 



—1 = [ (i|)-o) )cos0 + ^3 = (J> - zM 



C CWj. p 



C|4>ol c U, 



where M is given by 

M = S I i I + y WR^ i E ( 0 ) 

'o' S s TT O 

and ojp is the projectile angular velocity. 



(80) 



( 81 ) 
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Equations (66) , (67) and (80) are combined with Eq. (17) 

After simplification, and retaining only the first order 
(linear) terms, the results are: 



ilf = - (i|i,-cot0^) 0 



O' ~ 



<t> 

M ^O 



^ = (tJ^qCsc0q) 6 - c 



(82) 

(83) 



<t>. 



0 = 7^ 9 - i{JoSin0o [4> (1+X) -0Xi|i_sin0» + i{;Xcos0_] (84) 

- ° . o o . o 



These equations are then integrated under the following 
initial conditions 







= 0 


^(0) = 0 


> 


0 

II 

1 

0 


0(0) = 0 
o 




0(0) = 0 


0(0) =0 J 




0(0) = 0 


Eq. (82) yields 






^ = -(i|»QCot0^) 0 






From Eq. (83) 






i = -4>o (io^°®®o)® 


- “t 


^o 


c 1 





- X 



(85) 



( 86 ) 

By substituting the expressions obtained in Eqs. (85) and (86) 
into Eq. (84), the differential equation for 0 becomes 



.73 



A Uo 1 



( 87 ) 



or, in a more familiar form. 



0 2C ’ ’ 

■-- + e + e = D + 



Rt 



n n 



( 88 ) 



• • 

where % = '/'q = <l>o/l'i>ol = 



and 



Z = 



M 



2aUoI^o 



= sine = 

• o 



'Po 



- jSin 26 o 



= - li+AI &ine = - 






MD 

ioC 



Expansion of the expressions above, using Eq. (81) shows 
that the damping coefficient is given by: 



Z = 



ts|*J . U,WR3|E(e„)l 



or 



; = — J- [S + VsRsW I BOq)] 

2A<P„ ® ” Xif'oCosSo 



and the ramp slope (rate of decay) becomes 



7,4 







I 



5;. _ 




( 90 ) 






The integral of Eq. (87) with the initial conditions 
specified [0(0) = 0(0) = 0] is Eq. (65) where 

Js 



<I>T, = 2(j>^ = 2tan 
R D 



-1 (i-?M 



(91) 



For large values of ?Wj^t, the ramp decay dominates the 
response and the solution may be further approximated by 



0=0^ + 0=0^ + D + — (u t-2i;) 

Jn 



'O Z o ' n 



The time to nutate through a known angle 0 may be foxand 

from 



t 




2 ? _ 2 ? 



u 



n 



(0 



n 



(0 - 6 ) 
o 

R 




(92) 



where R/ D and are prescribed in Eqs. (88)-(90). 

Figure 15 is a plot of 0 vs. time comparing the results 
of the exact solution from Fig. 14 with those of the 
approximation developed here. The physical constants and 
initial conditions are identical for the two cases. 



D. THE EFFECT OF GEOMETRY ON OBTURATOR RESPONSE 

In the discussion to this point a primary consideration 
has been to model the external forces acting on the obturator. 
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Figure 15. Comparison plot of 0 vs. time for the (a) exact and (b) 
approximate solutions. 



The effect of geometry has not been discussed, but is a most 
significant contributor to the response of the obturator. 

C and A are functions of mass density and physical di- 
mensions. By varying the material composition (density) 
the response will vary. Utilizing the physical dimensions 
of the obturator in Table i, the obturator volume may be 
calculated by subtracting the volvimes of a right circular 
cylinder of radius r , length 2RgCOsSg and the two spherical 
caps at the ends of the hole from the volume of a sphere [26]* 
Thus the volume of the obturator is given by 

4 2 

V = -^R ^ (R -R cos8) ^ (2R + R cosB) ~ 2irr_^R^cos3 

Collecting terms gives 

V = 2irRgCos6 [Rg^ (1 - - rp^] 

By knowing the mass density, the weight of the obturator 
[for use in the torque expression, Eq, (61) ] may be obtained. 
Therefore as the density varies, the values of C, A and W 
vary. Fig. 16 is a family of curves of 6 vs. t for varying 
density with all other physical properties being those of 
Table 1. 

There is however a limiting case for the size of the 
hole. If the hole becomes too large (assuming same size 
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Figure 16. Comparison plot of 0 vs. time for obturator material 
densities of (a) 6000 Kg/m^ and (b)lOOOO Kg/m^ 



projectile hole) , there may not be sufficient obturator 
surface area to block the passage. The limiting relation- 
ship is 



r 

s 






R 



s 






This limit is purely theoretical and is not a feasible design 
parameter. Engineering judgement in the selection of realis- 
tic values is required so that the projectile will meet the 
ballistics criteria and retain its integrity while undergoing 
the forces associated with the entire firing evolution from 
chambering to target impact. 

The effect on the obturator response due to increasing the 
hole diameter is shown in Fig. 17 • The inertial contribution 
to the response (coupled term) is increased by increasing C/A 
and therefore X; however the friction torque (Mgg) contribution 
is lessened by the reduction in the mass of the ball. 
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Figure 17. Comparison plot of 0 vs. time for obturator hole radii of 
(a) 3.75 mm and (b) 5.58 mm. 



III. EXPERIMENTATION 



A. EXPERIMENTAL APPARATUS 

An experimental system was designed to simulate the 
spinning EOT. The apparatus consisted of a compressed- air- 
driven spin-up rig, optical timing mechanism, air manifold 
and associated piping, tubing and electronics. 

The obturator was fabricated from a standard 5/8- in. 
diameter chrome -steel (52100) bearing ball. The ball was 
annealed to allow machining and then bored along a diametrical 
axis (Fig. 18) . The bored ball was then mounted in a bakelite 
metalographic specimen mount. The mounted ball was placed in a 
milling machine and a flat was machined in the bakelite at a 
specified angle relative to the axis through the hole in the 
ball (Figs. 19 and 20) • This flat was then used as the 
polishing plane for metalographic specimen preparation. A 
small flat spot was polished on the obturator to provide a 
highly reflective surface at a known orientation relative to 
to the z-axis of the obturator (Fig. 21) . 

After removal from the bakelite mount, the obturator 
was placed in a three-piece, lucite housing (Fig. 22). The 
mating ends of the two hollow inner cylinders were each 
machined with a 5/8- inch end mill to a depth of approximately 
5/16-inch. When mated, a spherical cavity was formed to 
accomodate the obturator. The third cylinder was press fit 
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, Figure 18. 



Boring of hole along a diametrical axis of 
the annealed bearing ball. 




Figure 19. Machining of flat in metalographic specimen 
mount. 



82 








Figure 20. Machined specimen mounts with varying 
angles ( 6^) . 




Figure 21. Polished spot orientation on obturator and 
relationship to 9^. 
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Figure 22. Obturator and components of lucite obturator 
housing with mating ends of inner cylinders 
machined to house the obturator. 










i«ii 





















Figure 23. Assembled obturator housing with aluminum 
end pieces. 
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over the others after the obturator was inserted to insure 
alignment of the inner pieces and provide rigidity (Fig. 23). 

The ends of this lucite assembly were then press fit into 
aluminum end-pieces similar to those in Fig. 23 . These end- 
pieces served as the shaft for the bearings and one also served 
as the prime mover (bucket wheel) for the apparatus. The shaft 
rotated in two ball bearings mounted in aluminum pillow blocks 
aligned on a rigid pedestal. The prime mover was a bucket 
wheel machined from a solid alumimm disk (^ig. 24). 

The bucket wheel was driven by compressed air supplied 
from an installed system through an air filter to a manifold 
and then through one of two Model-10 Kendall pressure regu- 
lators. From the regulator, the air passed through a flexible 
tube to a tee, each leg of which supplied a nozzle. These 
nozzles were mounted opposite one another on the pillow block 
in such a way as to allow the air jet to impinge upon the bucket 
wheel to cause rotation (Fig. 25). The speed of rotation was 
sensed by a Bentley Nevada Proximitor, Model 3100N, which was 
mounted above the bucket wheel to detect the passage of each 
point on the wheel. The sensor was supplied from a LAMBDA 
Regulator Power Supply Model LP413FM by 18VDC. The pulses 
generated by the sensor were cotinted, averaged and displayed 
as a frequency by a Monsanto Programmable Counter-Timer Model 
HOB. 

While the spin-up rig was being brought up to the desired 
speed, an air jet from a nozzle mounted rigidly on the pedestal 



85 







86 



Figure 24. Assei^led spin-up rig with obturator oriented in 
loading" position. 
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Figure 25. Spin-up rig showing bucket-wheel with air je 
orientation and location of Proximitor above 
bucket wheel. 



at the end opposite the bucket wheel held the obturator 
fixed to the spinning rig. This jet was supplied from the 
common air manifold and the other pressure regulator also 
through flexible tubing. The air jet passed through a hole 
in the plunger of a Rocker Solenoid, R.S. No. 10-207. This 
hole was aligned with the hole through the "projectile" to 
allow the jet to impinge upon the obturator (Fig. 26) . When 
the switch at the far right of Fig. 26 was activated, standard 
115VAC was applied to the Rocker solenoid and removed from a 
normally closed ASCO Solenoid Valve in the air supply line to 
the jet. The plunger retracted from the position shown in 
Fig. 26 and the solenoid valve closed. The retraction of the 
plvmger performed three fxinctions. It first caused a pulse 
to be generated by another Proximitor movmted next to the 
nozzle and powered from the same 18VDC supply. This pulse 
started the timer function of a second Model IlOB Counter-Timer. 
The plunger also covered the nozzle outlet by misaligning its 
hole with that of the nozzle. This removed the restraining in- 
fluence of the air jet on the obturator and prevented any air 
remaining in the supply line from impinging on the obturator. 

And lastly, in the fully retracted position, the machined and 
polished end of the plunger was positioned opposite the pro- 
jectile hole (Fig. 27) . When the released obturator nutated 
through enough of an angle, the beam of a SPECTRA PHYSICS MODEL 
132 LASER, MODEL NO. 3187, passed through the projectile and 
was reflected by the polished end of the plunger into a light 
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Figure 26. Holding- jet nozzle (1) , Proximitor (2) and Solenoid 
plunger (3) extended prior to obturator release, 
(hole through plunger is aligned with the end of 








90 



sensitive diode (Fig. 28) . The pulse emitted by the diode 
triggered the stop channel cn the timer and the elapsed time 
was displayed. Figure 29 shows the entire experimental set-up 
with air flasks in the background. 

B. EXPERIMENTAL PROCEDURE 

The obturator was positioned in the spin-up rig in such 
a way that the laser beam was reflected off the polished 
spot back to the target mounted on the laser (Figs. 29 and 
30) . The center of the target is the location of the beam 
and the obturator was adjusted to place the reflected spot 
on the center of the target. Thus the orientation of the 
z-axis through the obturator was known relative to the Z-axis 
(laser beam) through the projectile. This is the initial 
value 0Q. 

The solenoid plunger was checked in position with the 
hole aligned with the nozzle (Fig. 26) and then air was 
admitted to hold the obturator in place. The pressure was 
then increased to the turbine to set the apparatus in motion. 
Figure 31 illustrates the circle traced by the reflected spot 
(Fig. 30) on the target. This circle facilitates alignment 
of the laser to insure the beam was centered through the pro- 
jectile. (It can be seen in Fig. 31 that the beam was slightly 
off-center.) Alignment was accomplished by positioning the 
beam from side to side for horizontal alignment, then raising 
or lowering the spin-up rig's mounting platform with adjusting 
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Figure 28. Light sensitive diode housing (1) illuminated by 
reflected laser beam from the polished end of th 
solenoid plunger. 
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Figure 29. Complete experimental apparatus layout. 





Figure 30. Laser reflection on the target from the 
polished spot on the obturator. 




Figure 31. Circle traced by the reflected laser while 
the apparatus is spinning. 
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screws for vertical alignment. When aligned, the circular 
trace of the reflection exactly follows the circles of the 
target. 

Once the alignment had been checked and the obturator 
positioned at the known 0^, the pressure was increased to 
the turbine nozzles until the desired spin rate was achieved. 
The regulator was then adjusted so as to maintain the spin 
rate at ±5 Hz on the digital display or approximately ±1.75 
rad/sec. (The actual tubular projectile spin rate is 
approximately 1.2xl0‘* rad/sec) . Once the desired rate was 
achieved and noted, and the timer checked and reset if neces- 
sary, the switch (Fig. 26) was activated. This secured the 
air to the nozzle holding the obturator and retracted the 
solenoid plunger (Fig. 27) . The obturator was released and 
began to move relative to the spin-up rig. When the angle 
9 reached the value at which the laser beam could pass through 
the obturator hole (0=9(,j-=8) , the beam struck the polished 
plunger and was reflected into the light sensitive diode 
housing (Fig. 28) . The reflected beam striking the diode 
caused the diode to emit a pulse which turned off the timer 
gate. The elapsed time from plunger retraction to release 
the obturator, to the obturator nutating to, 9^^ was dis- 
played on the timer display to the 0.0001 second. This time, 

t^^, was recorded, 
cr 
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C. DATA REDUCTION 



The only data reduction required was the conversion from 
Hz to Radians per second. The frequency read-out on the 
digital display indicated the number of teeth on the bucket- 
wheel passing the sensor in one second. Since there were 
eighteen teeth on the bucket wheel, the conversion was 
simply 
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Appendix B is a listing of the raw data and the 



corresponding values of Wp. 
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IV. CORRELATION OF RESULTS 



Figures 32 through 37 show 6 vs. time for various values 
of o)p from 40 to 1000 sec"^. At high and low values of Wp, 

0 approaches zero more slowly than at the intermediate values. 
Figure 38 is a plot of Eq. (92) and the experimental data. 

The plot of Eq. (92) uses the physical qualities from Table 1 
and y =0.35 over the range 100-^_-^1000 sec”^. It is clear 

S 

that the experimental data falls below the approximation curve 
at high Wp and rises sharply above the curve at low Wp. 

In order to determine what parameters may effect the 
position of the approximate curve, sensitivity checks were 
run on the quantities in which the most uncertainty existed, 
namely Ps and 6^ (even though 0q could only vary by ±1°) . 
Figures 39 and 40 indicate the sensitivity of t^^ to variations 
in Pg and 0^ respectively. Small variations in 0^ are not 
discernable within the range of scatter of the data (Fig. 41). 
Variations in Pg have a much more pronounced effect. From 
Fig. 42 it may be seen that the data at higher values of co 

ir 

correlate well for a Pg=0.55. 

It has been noted that the approximate solution is based 
upon the assumption M/CWp^<<l. At values of Wp(<200 sec~Mf 
this assumption is no longer valid and therefore calculations 
and correlation of the data to the approximate solution will 
not be conducted for Wp 200 sec~^. Also at low values of cOp, 
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Figure 32. 0 vs. time for gravity model, to = 40 sec 
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Figure 33. 6 vs. time for gravity model, w = 100 sec 
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Figure 34. 6 vs. time for gravity model, w = 200 sec 
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Figure 35. 0 vs. time for gravity model, (O- = 500 sec 
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Figure 36. 6 vs. time for gravity model, w - 700 sec 
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Figure 37. 0 vs. time for gravity model, w = 1000 sec 
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Figure 38. Experimental data and t^,j- vs. oip from Eq. (92) for 
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Figure 39. Plot of t^^ vs. Wp from Eq. (92) for values of 
Mg = 0.1; 0.25; 0.35; 0.55 
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Figure 40. Plot of t vs. o) for various values of 
^ or p 
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Figure 41. Plot of experiraental data and t^^. from Eq. (92) for 
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Figure 42. Experimental data and t^^ from Eq. (92) for 
u_ = 0.10, 0.25, 0.35 and 0.55 



t appears to approach «» almost assymptotically . Because 
cr 

00-90® (87.25®) , the obturator may be near a metastable state 
at low (Op such that the coupled terms in the equations of 
motion have little effect. The amplitude of oscillation at 
low (Op is also relatively small. 

In general, for large values of ?“j^t (Eq. 65) and 
neglecting the effect of viscosity, t becomes 



t 



X+1 M 

^ A(0p ^ COS 9 q 



Msin0Q 



[(00-6)- |sin6o] 



(93) 



where M is the gravitational contribution from Eq. (81) . 
This may be written as 



t = 




+ T (0 
2 



P 



(94) 



with T and T being constants determined from the known 
1 2 

physical parameters (Table 1) . Utilizing values of 

0Q = 87.25®, 0 = = sin"^ (r/Rg) = 36.7 34® and the measured 

weight of the obturator, W = 0.0798 N, 

M is found to be 

M = 6.34x10”'* y N-m 



There fore , 

T^ = 2.039x10^ y 



and 



T = 3.03x10"“ / y (95) 

2 
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thus 



t ^ 2.039X1QS - ^ 



3.03x10 



-4 



CDr 






But for large values of 



3.03xl0-‘*(o 
t = P 

V 

A least squares fit of the data Cc0p^35.97) to Eq. 

such that t = T cOq, gives 

2 ^ 

2t, 2 

T = yT 

2 2tiC0pi 

Substituting in the data points (denoted by ' * ' in 

Et^^ = 2.1046884 sec^ 

Ewpi^ = 2.8413714xl0®sec“^ 

Eti^w . 2 = 2.430257x10^ 

^ pi 

there fore 

T = 8.6603508x10““ sec^ 

Since the standard deviation for T is given by 



= 1 



Et,. 



(j2.p4 

2 
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(96) 

C96) / 

Appendix B ) , 




m 

i 





where 



2 1 
^ = N 






pi 






the following results are obtained: 



= 1.4659547x10^ 

= 3.9181053x10"^“ 

a„ = 1.97742x10"= 

T 



^ =0.023 

T 

2 

The data therefore, deviates from the least squares value 
by only 2.3 percent. 

Using the value calculated for T in Eq. (95) , 

2 



3.03x10 
8.66x10"** 



with EC87.25°) = 1.00458, 




1 (1.00448) = 0.63954 

TT 



thus 



0.35 

Us = D"’639B4 = 0.5473 
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with this value of y , values of t were obtained from 

s cr 

the exact solution over the range 125^a>p^700 sec“^ and 
plotted with the approximation of t^^. from Eq. 94 and the data 
points CFig. 43). 
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V. CONCLUSIONS 



The response of the obturator (9) as a function of time 
to applied torques was found to be that of a rigid body of 
revolution about its mass center. Euler's Modified Equations 
of Motion in terms of Euler Angles describe the response to 
these applied torques. The response is predictable and was 
compared for various combinations of physical parameters. By 
varying density, hole size and angular velocity, various re- 
sponse times may be obtained. It was found that in each 
case, a minimum response time was attainable for a given value 
of Up. For the case of the experimental apparatus, this mini- 
mum occured at u = 150 sec“^. 

P 

The analytical model and a linear approximation to this 
model were used with the parameters of the experimental 
apparatus and close agreement to the data was achieved at 
u^>200 sec“^. Below Up = 200 sec“^ the exact solution appeared 
to be higher and the approximate solution lower than the experi- 
mentally observed values of t^^j-. This is explainable in the 
approximation since it was based upon the assumption that 
M/Up^C<<l. At low values of Up, this assumption does not hold. 

The trend in the exact solution at low Up is to approach 
tpp = “ almost assymptotically . This trend is observable in 
the experimental data as well and may be explained by the 
metastable initial position of the obturator 9p-90° (6o=90° is 
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a stable position). The effect of the coupled terms in the 
equations of motion is small at low therefore the initial 

effects are not as great and the amplitude of oscillation of 
0 is insufficient to cause departure from this metastable 
position. Since the actual projectile operates at o)p=12,000 
sec”^, high values of o)p are of interest. The theory agrees 
well with the data above o)p=200 sec“* and therefore the 
correlation to the data was made at these higher values. 

Because the linear approximation agreed closely at higher 
u)p, a least squares fit of that approximate equation to the 
experimental data (o)p— 235.97 sec“M lead to the calculation 
of a higher sliding friction coefficient than had been used 
previously. The assumed value was iig=0.35 for laminated 
plastic on steel [27]. This was found to be too low by the 
least squares fit. A sensitivity check of the approximate 
solution for t^^j- to variations in Pg also indicated that 0.35 
was too low (Fig. 39) . Therefore, a method for determining 
the sliding friction coefficient from the experimental data 
and the approximate solution was discovered. 

By knowing initial orientation (Sq) / projectile spin rate 
(ojp) and the physical dimensions and material properties of 
the projectile, the response may be predicted in a gravity 
environment with reasonable accuracy. 

The aerodynamic model appears to give reasonable results 
within the scope of this study. The complex forces acting on 
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the obturator in flight are extremely difficult to model 
theoretically and therefore they were simplified greatly. 
The confidence gained in the gravity model, however, leads 
to the conclusion that once the complex aerodynamic forces 
are better understood, the motion of the ball within the 
flight model will be predictable. 
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VI . RECOMMENDATIONS 



Further experimentation should be undertaken with 
variation of other parameters (hole size, obturator material) 
to obtain more confidence in the analytical model. Wind 
tunnel tests of the EOT should be conducted in order to more 
closely approximate the actual forces on the obturator in 
flight, and to determine the position (0) at which the de- 
tached bow shock is swallowed. This would greatly enhance 
the predictive capabilities of the model and thus allow 
coupling to a trajectory model for predicting flight char- 
acteristics. This, in turn, could be coupled to a numerical 
optimization routine for optimizing the EOT design for the 
desired performance criteria. 
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APPENDIX A - COMPUTER CODE 



4(4e 9C itx :di4i v 4i A4i *« 4c*c* 

BALL OBTUPAT'-D TUBULAR PRQJf-CTILn » 



BALL MOTION ANALYSIS 
LT J* W* BLOOMER lit USN 



¥ 



|i ^I|ii«c 4 c 4 « * « « :* « A* Ik 

DIMENSION XOO) tXODTOQI ,C( 15) 

C(l)s (SI N(RCT/ ) )=>^-‘2t WR'-R*’ Bl^TA IS GIV^N BY ARSIN(RP/RO> 

C(2>- LAM8CA = (C-A)/A 

C(3)= a;kCOYN;MIC np GRAVITATIONAL SLIDING FRICTION CONSTANT. 
C(4)= total KINrTIC ENERGY OF CDMBINFD PROJlCTILr ANO OBTURATOR. 
C(5)= MIM*^LM Off !N*^ctIA» A. 

CC6)= MAXLVUM PCM-NT OF U'5RTIAt C. 

C(7)a MCKGNT CF IM^ffTTA C »= A FOLLOW CIRCULAR CYLINDER 
approximating that Zf THE PROJCCTILE# 

C(8)* FLUID FRICTICN CGNSTAMTt S. 

CC 1U)=1. 

I 



c 

8 
c 
c 
c 
c 
c 

c 

8 
c 
c 
c 
c 
c 
c 

c 
c 

C*««^»##CALL INTEG2 SIMULTANFCUS OIFF. EON. SOLVER 
C 

c 

I CALL INTcG2(TtX,X00T,C ) 

C 

C 

ELLER ANGLES ThETAt PHI ANO PSI AND 1ST OERIVATIVTS 
C 
C 

TH=X(1) 

X0CT(1)=XC2) 

PS=X(3) 

X0CT(3 )=XC4) 

PH»X(5) 

XCCT(5)=X(6) 

C 

C ♦ 

TRIGCNCMETRIC FUNCTIONS OF EULFR ANGLES 
C 

c 

STh=SlN(TH) 

CTF=COS(Th) 

SPH=SIN(FH) 

CPF=COS(FF ) 

SPS^SINC PS ) 

CPS=CJS(PS) 

C 

TFETA-DQTt PHI-CCT AND PSI-DOT 
C 

THC^XDGTd) 

PSD=XOGT(3) 

PHO*XDOTC5) 



^^*^^***CALCULAT5 ANGULAR V'"'LGCITIES WITH R-SPECT TO BODY- 
C FIXED AXcS At 8t ANO C RE SPCCTI VCL v. 

C 

□ Ml-PSD ’'STh + SFH^THD^^C PH 
aM2-PSD*STh<CFH-THD^SPH 
CM3=PS O^CTh^FFC 
C 

C***^^^=^CALCULATi ANGULAR VELOCITIES ABOUT INERTIAL AXCS X, Yt 
C 2 RcSPlCTUGLY. 

c ■ 

OMXsGMl* (CPH*rPS-SPH^CTH>*£PS )-fCW2*(-SPH*CPS-CPH*CIH^SPS) 

♦ ♦CM3^(STF>'5FS ) 

OMY^OMl-^ (CPH^SPS«^SPH«CTH^CPS KCM2* (-SPH« S PS + CPH*CTH^CPS ) 

♦ 40K3’^(-STh + CPS ) 

QMZ=CMl*STh^SPH^CM2*STH^CPH4CM3«CTH 

S*=*»***»*CALCUL4TF Th’^ 4NRULAR VTLCCITY OF TH:r PROJ5CTIL& 

C FROM CCNSiFVATION OF -N.'RGY PRIMCIPLcS. 

C 



AMO 










m 







SPINsO.O 



SP1N*S0<^T( (2.*C(4J-C<5)*(CMl**2+aM2**2)-C(6)»OM3*«2)/C(7) ) 

C 

C*-###ii*****CALCULATc angular vlocity of ball relative to thf proj^ctil 

C WITH PESP-:CT TO BOCY-PIXEO AX^S. 

c 

OMIR »aMl-SPlN*SPH*^TH 
QP2RaCP2-Sf IN»CPI-*'>TH 
OM3R»Of'3-SFIN*CTH 
C 

C***#****CALCULATF ficLATIV* ANGULAR V-LOCITI'S WITH RESPECT TO 
C THc INIRTIAL P-F! RCNCr FRAM^. 

^ 0HXR30f‘lP=»<CFH'‘CPS-SPH*CTH’»SPS I ♦0M2R*(-SPH*CPS-CPH»CTH*SPS) 
•+OH3R*(STH*SPS ) 

Q»‘YRsOM1R»(CPH*5PS+SPH*CTH*CPS)*OM2R*<-SPH*SPS+CPH*CTH«C'>SI 

♦♦CH3R*<-STH«CFS) 

0MZKs0mR*STh*S?H+nM2R^STH»CPH+0M3R*CTH 



C 

C***s,^»*:^CALCULATE HAGNITU02: CF RELATIV*! ANGULAR VELOCITY VECTOR. 

C ' 

OPR=SQRT (CN1F*Cmir+0H2R*CH2R4CM3R*0M3R) 

C 

C*#**»***CALCULATE TH= MOMENT aPM AT WHICH THS FRICTIONAL FORCFS ACT. 

C 

C ARMa MOMENT ARM AND C(l» tS AS OCFINFD EXCEPT IN THE CASS OF 
C GRAVITY WH«=RC C(H = 1.). 

C FRICa CGMBINcC FLUID AMD SLIDING FRICTION CCNSTANTS. 

C 

ARPaSORT (l.-C (1 >«(OMXR/CMR)**2-U.-C (1) J *(QM2R/CMR)**2) 
FRICaC(3»aARM 
AMC.M=C(8 >aCHP4FRIC 
C 

C***»J>:a**#aCALCULATr MCMHNTS ACTING ON THE BALL ABOUT TH= BODY-FIXED 
C AXES. 

C 

Gla-AM0M*GM1R/CMR 
, G2*-AMCH*CN2F/CMR 
*'• G3a-AMCH»CY3?/C“R 

Ca**»*****aCALCULATS .*LSM;NTS OF MATRIX OF COEFFICIENTS. Yl, Y2 ANC Y3 . 
C 

YlaGl/C( 5)-FS0«r^H0PCTH*SPH4PHC-'STH*CPH) ♦TH0*PHC*SPH-C(2)» 

♦ (PSD»STh4CPH-THD4SfM)*( PSD'-C 7H+PHO) 

Y2=G2/C(5)-FSC-» ("HC“CTi-«CDH-PHO*STH*SPH)4PHC*THC*CPH4C( 2) * 

*( PH04PSP*C1H »a(PSD •SrHa*:PH+THO-*CPH» 

Y3=G3/C(6) «FSCa'>0*5TH 
C 

C**a***p*a,^CALCULAT.S ANGULAR ACCHLSRATION TFRMS FOR INTEGRATION, 2ND 
C Of RIVATIVFS OF "ULTR ANGL'S WRT TIME. 

C 

XOOT(2 » = CPI-'Y1-S?H’'V2 

XCCT(6 > = -SFH*CTH/STHaYl-C^H*CPH/STH*Y24Y3 
XDCT(4)=SPh/STH«Yl4CPH/STHaY2 
C 
C 

C»»*a^»*CALCULATS VALUES FOR THE APPROXIMATE SOLUTION 
C 

c 
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THCRITa.6411 

IF(T.N'c.C.CCCCOCO) GO TO 25 
THNCT=X(1 ) 

OMPaXC 4) 

ERTCaC(6 ) 

AaC(5» 

AHCa=C (2) 

AHCAP=AHCA41 . . 

ENaC (3 )>i'-. 63S5A 
Da-AMDA/2.-*SIN(2.4THN0T) 
Ra-AMCAP^iNaSINCTHMOT)/ ( -P'^C^OMP) 

TCRIT=-1 ./R4( THNOT-ThCRI' )-D/R 
WRITF(6,2>> TCRTT 
FCPMAn5X,‘fCR!T = '.SIB.B) 

ZKTA=AMDAP/(2. •‘A^^AMDAaCN PaCMPaCCS(THNOT) MCH 
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I 
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IFUHTA.Gr.l^OOn TO 25 
OZtETA=SCHT (1 

PHAS^=ATAN(CZ':TA/CZ-TA^ (1.-2 ./ (A ZHT Ai«Z ^'T A-AMDA P ) ) ) ) 
THPR0X=ThNGT+C+P/aMP'^(-2.^ZtTA)+l./Z^TA*SIN(PHASIf) 

GO TC 26 

25 THP.-<OX=THNr:T>0>R/GMP«(CMP^T-2.'*Z[:TA>>l./ (EXP ( 2€T A^OMP^T )^DZ r'TA I ♦ 
<«SIN(OMP=^T»CZc:TA>PHAS<:) 

26 X(7)=THPRCX 
XC9>=OMZ 
XU3 )=SPIN 
GO TO 1 
END 
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APPENDIX B - RAW DATA 



READOUT 


TIME 


“p 


READOUT 


Tir-iE 


“P 


[HzJ 


I secj 


[sec“^ J 


[HzJ 


[see] 


[sec“ ^ ] 


776 


. 274 


270.88* 


880 


.2717 


307.2* 


6 76 


.218 


235.97* 


920 


.3026 


321.14* 


679 


.2181 


237.02* 


920 


.2847 


321.14* 


573 


.2157 


200.01 


936 


.2816 


326.7* 


570 


.2079 


198.97 


983 


.2699 


343.1* 


574 


.2009 


200.36 


1001 


.306 


349.4* 


570 


.19 7 


198.97 


994 


.2972 


346.97* 


571 


.2031 


199.32 


1055 


.296 


368.26* 


475 


.2099 


165.8 


1063 


.2530 


371.06* 


479 


.1833 


16 7.2 


1151 


.315 


401. 77* 


570 


.2179 


198.97 


1022 


.369 


356.7 * 


475 


.215 


165.8 


1175 


.378 


410. 15* 


479 


.223 


16 7.2 


1258 


. 3348 


439. 12* 


470 


.2 30 


164. 1 


1278 


.3968 


446.1 * 


428 


.1831 


149.4 


1255 


.4268 


438.08* 


425 


.1786 


148.35 


1221 


. 3088 


426.2 * 


430 


.1659 


150.1 


1170 


.284 


408.4 * 


367 


.1896 


128.1 


1158 


.281 


404.2 * 


363 


.2037 


126.7 








370 


.19 59 


129.15 








277 


.2249 


96.7 








267 


.2292 


93.2 








279 


.2259 


97.4 








323 


.2 30 


112.7 








323 • 


.2149 


112.7 








225 


.493 


78.5 








223 


.658 


77.8 








223 


.599 


77.8 








226 


. 3241 


78.9 








807 


.2592 


281.7* 








810 


.2339 


282.7* 








853 


.2773 


297.8* 








853 


.2631 


297.8* 








850 


.2917 


296.7* 








854 


.267 


298.1* 








921 


. 3022 


321.49* 
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